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Abstract. We extend the work of Kurchan on the Gallavotti-Cohen fluctuation 
theorem, which yields a symmetry property of the large deviation function, to 
general Markov processes. These include jump processes describing the evolution 
of stochastic lattice gases driven in the bulk or through particle reservoirs, gen- 
eral diffusive processes in physical and/or velocity space, as well as Hamiltonian 
systems with stochastic boundary conditions. For dynamics satisfying local de- 
tailed balance we establish a link between the action functional of the fluctuation 
theorem and the entropy production. This gives, in the linear regime, an alterna- 
tive derivation of the Green-Kubo formula and the Onsager reciprocity relations. 
In the nonlinear regime consequences of the new symmetry are harder to come 
by and the large derivation functional difficult to compute. For the asymmetric 
simple exclusion process the latter is determined explicitly using the Bethe ansatz 
in the limit of large N. 



1 Introduction 



The study of stationary nonequilibrium states (SNS) of macroscopic systems 
evolving according to classical (or even quantum) mechanics, which are kept 
out of equilibrium by contact with thermal reservoirs has a long history. The 
interactions of the system with the reservoirs are usually modeled by the addition 
of stochastic boundary terms to the deterministic evolution equations describing 
the isolated system [|]]-@. Fully deterministic evolutions of systems coupled to 
infinitely extended reservoirs have also been investigated ||-[10|, as have fully 



stochastic models [jnj-fllj]. A more recent development has been the study of SNS 
of systems evolving via deterministic "thermostatted" dynamics [1^, |16|, [IT]. In 
this scheme the SNS is maintained by deterministic external driving forces which 
are balanced by a "friction term" suitably chosen so that the system evolves on a 
compact surface (generally one of constant energy) in phase space. The resulting 
dynamics no longer conserves phase space volume, which is contracted on the 
average fl9f . Consequently, the SNS is singular with respect to the Liouville 
(volume) measure induced on the surface. 

While the jury is still out on what has or can be achieved by such a dynam- 
ical systems approach, there are some intriguing theoretical results which have 
emerged |2(], ^TJ . In particular Gallavotti and Cohen [^] motivated by results of 
computer simulations [^3j discovered that under suitable assumptions these SNS 
satisfy a certain symmetry, which they doobed a "fluctuation theorem" . Assum- 
ing that the dynamics satisfies time reversal invariance and is sufficiently chaotic, 
so that the SNS is given by an SRB measure, they prove that the probability dis- 
tribution for the phase space contraction averaged along a trajectory over the 
time span r has, for large r, a highly non-obvious symmetry, whose specific form 
will be given below. Near equilibrium the fluctuation theorem implies Onsager 
reciprocity and the Einstein relations. We refer to ]24j for an exhaustive discus- 
sion of the theorem and of the "chaotic hypothesis", which implies the validity 
of the theorem for more general SNS. 

In a recent work Kurchan showed that, with proper definitions, the fluctuation 



theorem is valid also for certain diffusion processes [25]. In the present work we 



extend Kurchan's results and show that for Markov processes the fluctuation 
theorem holds in great generality. Our proof is based on the Perron-Frobenius 



theorem and goes much beyond the cases considered in [25|. While the assumed 



stochasticity of the dynamics should provide "for free" a suitable modification of 
the chaotic hypothesis it is far from obvious how to find the quantity that plays 
the role of the phase space contraction. In analogy to the thermostatted case 
it is reasonable to expect that it will have some relation to the production of 
Gibbs entropy. In fact, this comes out from our general construction under the 
additional assumption of local detailed balance. This condition seems to play, in 
stochastic dynamics, the same role as time reversal invariance in deterministic 
thermostatted dynamics. When local detailed balance is satisfied, then the action 
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functional is proportional to the current of the conserved quantity, e.g. in the 
case of bulk driven lattice gases it is the bulk particle current. If local detailed 
balance is not satisfied, then the action functional can still be defined but it does 
not allow for such a direct physical interpretation. 

Near equilibrium the fluctuation theorem, generalized to systems with several 
currents, yields the Onsager symmetry and the usual Kubo formulae for the linear 
transport coefficients. These give a relation between linear response and current 
fluctuations in equilibrium. In this sense the fluctuation theorem, which holds 
also far from equilibrium, can be thought of as a natural extension, in the spirit 
of GC, of the fluctuation-dissipation theorem which holds for systems close to 
(local) equilibrium. In fact, when the large deviation functional is approximately 
quadratic over a sufficiently large range of the driving field, the fluctuation the- 
orem does indeed yield such an approximate extension. We do not, however, 
have a quantitative criteria of when this approximation is valid. Even aside from 
this extension it would be of great interest to have an experimentally verifiable 
consequence of the novel symmetry predicted by the fluctuation theorem for the 
distribution of large deviations in current carrying systems. This seems difficult 
to obtain for macroscopic systems since the time required for observing a large 
deviation in such a system far exceeds any observation time. 

The physical meaning of the fluctuation theorem can be understood best 
through the application to simple models. We therefore determine here the ac- 
tion functional for which such a fluctuation theorem holds in various examples: 
bulk and/or boundary driven lattice gases, the Fokker-Planck equation with driv- 
ing forces and a spatially varying temperature, and Hamiltonian particle systems 
driven by boundary reservoirs. In general, the rate function appearing in the fluc- 
tuation theorem cannot be computed explicitly and one has to rely on numerical 
simulations [|23| , [26|1 . In stochastic models the task is easier as long as we stick to 
noninteracting particles. We discuss one example of interacting particles, where 
the rate function can be computed via the Bethe ansatz. When the size of the 
system goes to infinity the rate function becomes singular. This is analogous to 
the behavior of thermodynamic potentials at a phase transition. This is done in 
the Appendix which may be read independently of the rest of the paper. 

Our analysis extends directly to discrete time stochastic processes, i.e. to 
probabilistic cellular automata. In this situation it is natural to think of the GC 
fluctuation theorem as a symmetry of discrete space-time Gibbs measures, as has 



been pointed out by C. Maes |[27|| . In fact, the framework of space-time Gibbs 
measures contains as special cases both the stochastic dynamics studied here and 
the thermostatted systems satisfying the chaotic hypothesis. In the latter case 
the space-time Gibbs measure is constructed through a Markov partition of phase 
space. 
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2 The fluctuation theorem for jump processes 



We start with the fluctuation theorem for a stochastic time evolution governed by 
a master equation on a state space with a finite number of points. This minimizes 
technical complications and at the same time provides a blue-print for stochastic 
systems with Langevin type dynamics. In this section we investigate the abstract 
structure, specific examples and applications will be discussed in later sections. 

2.1 Action functional, large deviations 

We consider a continuous time Markov jump process with finite state space S. 
Points in S are denoted by a. The jump process is determined by the rates, 
k{&, cr') > 0, for jumping from a to a'. More precisely, if the system is in the 
state a it waits a random time t > distributed according to the exponential 
law r(o")e~ r<CT )*<it, r(cr) = J2 a t k(a, a'), and then jumps to a' with probability 
r(cr) _1 A;(a, cr'), etc.. The generator acting on functions / : S — > K. is given by 

Lf{a)=Y,k^a')[f{a')-f{a)]. (2.1) 

cr' 

This means, if fi(cr,t) is the probability distribution of a at time t, then the rate 
of change of the average of /, (f)^(t) = E CT M°"> *)/(<r)» is S iven b Y 

j t (f),(t) = (Lf)» {t) . (2.2) 
Taking / to be a Kronecker delta at a (2.2) corresponds to the master equation 

= £ AV, t) - r(<7)M<7, = L*v(a, t) (2.3) 

cr' 

for the evolution of n(a,t), where L* denotes the adjoint of L. 

We assume that if k(a, a') > 0, then also k(a' ', cr) > 0, and that from every 
a all other a' G S can be reached by a succession of steps with nonzero rates. 
It follows then from the general theory of finite state Markov processes that 
there exists a unique stationary measure, /i s , which satisfies L*/i s (a) = and is 
strictly positive, // s (er) > 0. Starting from any initial state ctq, ji s is approached 
exponentially fast in time. 

We could as well take the state space S to be countable or replace S by W 1 . In 
the latter case one prescribes the rates k(x, x')dx' for a jump from x to the volume 
element x' + dx'. Setting r(x) = f k(x, x')dx', the waiting time at x is exponential 
with rate r(x) and the jump probability to x' + dx' is r(a;) _1 A;(x, x')dx' . Formally 
the difference is only in notation. However, now one has to make sure that the 
jump process has a unique invariant measure. 
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Let s — > a s be a trajectory (or history) of the jump process and define for it 
the quantity 

W(t,{a s ,0<s<t})= [ Yw a , a ,(s)ds (2.4) 



with w a y(s) a sequence of 5-functions, located exactly at those times s when a s 
jumps from a to a', with weight 

w(a, a') = log k(a, a') - log k(a', a) . (2.5) 

This means, if the trajectory o~ s , < s < t, visits in succession the states 
(To, <7i, a n , where ctq is the state at time and a n the one at time t, then 



W(t,{a s ,0< s<t}) =log 



fc(<J ,(7i) k(a n -i,a r , 



fc((Ji,(7o) fc((J n ,(7 n _i) 



(2.6) 



(Note that at an allowed transition /c(cr, o"') > and, by assumption, also 
£;((/, cr) > 0). 

For lack of a better name we call W(t), and the similar quantities to be 
defined below, an action functional. (We will generally not indicate explicitly the 
dependence of W(i) on the stochastic trajectory {a s ,0 < s < t}). 

Let (•) denote the expectation over all trajectories in the stationary process, 
i.e. starting in the steady state /i s (o"). We consider the generating function 
(exp[-AW(t)]) o f W(t) and define 

lim -ilog(e- w(t) ) = e(A). (2.7) 

t^oo t 

The fluctuation theorem we shall prove states that the limit (2.7) exists with e(A) 
convex downwards, and that 

e(A) = e(l-A). (2.8) 

By a general result, analogous to the equivalence of ensembles in equilibrium 
statistical mechanics, (2.7) implies a large deviation property for the probability 
distribution pt(w) of W(t)/t, i.e. for large t 

Pt (w)^e-^ w \ (2.9) 

where e is the Legendre transform of e. e{w) is convex up, e{w) > 0, e(wo) = 
for w = limt^oo^W^)), and by (2.8) e satisfies 

e(w) - e{-w) = -w , (2.10) 

i.e. the odd part of e is linear with slope —1/2. 
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To be more precise, from (2.7) one concludes that for every interval / 

lim --logProb(-W / (t) G I) = mine(w) (2.11) 

t^oo t t w£l 

with 

e(w) = max{e(A) — Xw} 

A 

= max{e(l — A) — Xw} = max{e(A) — (1 — X)w} 

A A 

= e{-w)-w. (2.12) 

Eq. (2.10) is of the same form as the GC fluctuation theorem, with the phase 
space contraction integrated along an orbit of the thermostatted dynamics in GC 
replaced by (2.5) summed over the jumps along the stochastic trajectory. We 
note that the convexity of e(A) and (2.8) imply that (w ) = (de(X)/dX)\\ =0 > 0. 
This is analogous to the result [|Tj| that the mean phase space volume contraction 
in the stationary state is nonnegative. 

To prove (2.8) we first define 

g(a,t)=E a [e- xw ^} (2.13) 

as the expectation value of e~ xw ^ conditioned on the system being in state a at 
time t = 0. We then have 

f t g(a,t) = ^A;(^cr')e- A ^^Kt)-r(^(a,t) 

a 1 

= v') 1 -^', a) x g(a', t) - r(a)g(a, t) 

= L x g(a,t). (2.14) 

Eq. (2.14) is to be solved subject to the initial condition g(a, 0) = 1. Therefore 
{e -xw it)) = J2fi s (a)g(a,t) = (e^V ■ ( 2 - 15 ) 

a a, a' 

Clearly (e LA ^) crcr / > by our assumption on k(a, a') above. The Perron- 
Frobenius theorem states that L\ has a unique maximal eigenvector fx which is 
characterized by f\(cr) > and satisfies L\f\(a) = —e(X)f\(a) with real maximal 
— e(A). This implies the existence of the limit in (2.7). Furthermore L\ and L* x 
have the same maximal eigenvalue. Using now the definition (2.16) we see that 

L*x = Li_x ■ (2.16) 

Hence the maximal eigenvector fx of L* x satisfies 

L*xfx(cr) = -e(X)f x (a) = In-JxW) ■ (2-17) 
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Since Li_a/i-a(c) = — e(l — A)/i_a(ct) and since /a > 0, we conclude by unique- 
ness that fx = fi-\ and e(A) = e(l — A). 

2.2 Time Reversal 

For a Markov chain with transition probability p(a, a') from a to a' the history 
{a} = {a , . . . , cr n } has the probability 

P({cr}) = fjL 8 [a )p(<To, «7i)...p(<7 n _i, <7 n ) (2.18) 

when starting in the stationary measure In analogy to (2.6) we define the 
action functional as 

n 

W(n,{a}) = -Y i log\p{a j ,a j - 1 )/p{a j - U a j )]. (2.19) 

i=i 

If we now denote the time reversed history of {a} by R{o~} = {cr n , ...,a }, then 
(2.19) can be rewritten as 

W(n,{a}) = -\og[P{R{a})/P{{a})]+]og\ Ma {a n )/fi 8 (a )] 

= W(n,{a})+log[pi s {a n )/fM s (a Q )). (2.20) 

One now immediately verifies that (e~ xw(jl ^) = ^e~^~ x ^ w ^) , where the average 
is with respect to P({<r}). Thus for W(n) the fluctuation theorem holds even for 
finite n. 

The form (2.20) is analogous to the form which the GC functional takes for 
thermostatted systems when the latter deterministic dynamical evolution is de- 
scribed via Markov partitions and it shows explicitly the role of time reversal in 
the original GC theorem. This observation forms the basis of the analysis of C. 



Maes in [27] , where he studies the fluctuation theorem for space-time Gibbs mea- 
sures and also discusses symmetry transformations different from time reversal. 

For the time continuous jump process a path {<j s , < s < t} is time reversed 
as {er t _ s ,0 < s < t}. The time reversed process has the jump rates k R (cr,cr f ) = 
k(a', a)fj, s (a') / /i s (cr) whereas the inverse waiting times, r(a), and the steady state, 
/x s (cr), remain unmodified. If P[ 0i t] denotes the path measure of the stationary 
process in the time window [0, t] and PR ^ the one of the corresponding time 
reversed process, then PE^ has a density relative to P[ j] and 

P[o,t] = e- w ^^ s {a t )/ii s {a Q ))P m (2.21) 

with W(t) from (2.4). Thus, up to boundary terms, the action functional equals 
— log(<iP[Q t ]/dP[o,t]) with dP R /dP denoting the Radon-Nikodym derivative. 

(2.21) remains meaningful for a stochastic dynamics which is not Markov and 
it can thus be used as a definition of the action functional in such a more general 
context. 
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One may wonder whether time-reversal could be replaced by some other 
transformation on path space. An obvious candidate is an internal symmetry 
S : S — > S such that S o S — 1. Using this transformation in (2.21) the GC 
action functional becomes 

W{t) = - f ds\og[r(Sa s )/r(a s )} (2.22) 
Jo 

up to boundary terms. From the examples we have studied it seems that only 
time-reversal leads to an action functional which has a simple physical interpre- 
tation. 

2.3 Sum of several generators 

We can generalize the above analysis to the case where the generator in (2.1) is 
the sum of several generators: L = X^li ^ 'with 

L {j) f(<j) = 5> W)[/(O - /(*)] . (2.23) 

a' 

We require that for transitions between a given pair a, a' at most one of the rates 
k^ is different from zero. It follows then from our assumptions on k(a, a') that 
if k^\a,a') > 0, then also k^\a',a) > 0. With this decomposition we define, in 
complete analogy with (2.4)-(2.6), the j-th action functional 

W^(t,{a s ,0 < s<t})= f (,)d. (2.24) 

(7,(7 

and the logarithm of the generating function 

1 m 

e(Ai, . . . , A m ) = lim -- log(exp[- V \jW^\t)]) . (2.25) 

i=i 

Repeating the arguments of Section 2.1 we conclude that the fluctuation theorem 
now takes the form 

e(A 1 ,...,A m ) =e(l-A 1 ,...,l-A ro ). (2.26) 
The corresponding analog of (2.10) is 

m 

e(twi, ...,w n )- e(-w u -w m ) = - wj (2.27) 

We shall use the general form (2.26) in Section 8 to prove the Onsager relations 
but note here that there are many ways of splitting L into a sum of L^'s. In fact 
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we can choose a different for every pair of points (er, a') for which k(a, a') ^ 
0. In all cases setting Xj = A for all j we recover (2.8) and (2.10) for W = 
^ r ^ =l W^\ We do not expect in general that the distribution of {l/t)W^\t) 
separately satisfies any fluctuation theorem. On the other hand we can consider, 
as a degenerate case of (2.23), a situation in which our system is composed of 
several independent systems with generators In that case, of course, each 
L^) satisfies the symmetry relation, ej(Xj) = ej(l — Xj), e.j(wj) —ej(—Wj) = —Wj, 
and e(Ai,...A m ) = EJLi e i( A i)> e(w!,...,w m ) = EJli^K). 

2.4 Entropy production 

For the probability distribution /i(a, t) at time t the Gibbs entropy is, as usual, 
given by 

S G (pt(t)) = -J^(<7,f) log n(a, t). (2.28) 
Using (2.3) the rate of change of Sq can be written in the form 

poW)) = R(n(t)) - A(n(t)) , (2.29) 

where 

RW)) - \ Y} k ^ OM^, *) - t)} bg[ ^ f^\ a '\ ] > o (2.30) 

Z aa , f J '\ <7 ' l ) K \U ) a ) 

and 

A(Kt)) = (I)Kth ^) = E*(^Olog[fe^]. (2.31) 



(When the generator is a sum, L = Y^j=i^\ as m ^ ne previous section, then 
the rate of change in Gibbs entropy is also a sum over j and consequently R = 

The breakup in (2.30) has the property that R is non-negative, while A is 
linear in /i. This suggests that we identify R as the entropy produced by the 
stochastic jumps and A with the entropy flow which can have either sign. Defining 
now the "integrand" in A, 1(c), as a "microscopic" entropy flux we note that the 
expectation value of W(t), starting with some initial distribution fj,(a, 0), is given 
by 

(W(t)) m = j ds(I) Ks) . (2.32) 
Jo 

We can therefore say that the average action equals the entropy flow integrated 
over the time span t. There is of course some arbitrariness in this identification. 
We could for example add some positive linear term to R and subtract it from A. 
Our choice of splitting is however quite natural in the examples discussed below. 
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When the system is in the steady state fi s , ^Sg{Hs) = and the entropy flow 
balances the entropy production, i.e. 

R(pt s ) = A(pt s ). (2.33) 

Thus 

{r l W{t)) = R(n a ) = A{n s ) (2.34) 

and hence the average action equals the rate of entropy production in the steady 
state. This suggests to identify W with the microscopic production of Gibbs en- 
tropy, a role played by the phase space contraction in the context of thermostatted 
systems; see also Section 7. 

We remark that R(/i s ) = A(/i s ) = in case k(cr, a') satisfies detailed balance 
with respect to the invariant measure /j> s (ct) ~ e~ v( - a \ i.e. 

e- y ^k(a, a') = e^^k^a' , a) . (2.35) 

Then each term in the sum defining R(fi s ) in (2.30) vanishes. This implies by 
(2.29) that A(fi s ) = (7) Ms = 0. The action functional (2.6) now becomes a 
telescoping sum with 

W(t) = V(a ) - V{a t ) . (2.36) 

It then follows from (2.7) that e(A) = and the fluctuation theorem becomes 
empty. More generally if 

k(a, a') = k (a, a')h(a, a') (2.37) 

and ko(cr, a') satisfies detailed balance while k\ is a nonequilibrium driving "force" , 
then the log/c terms of W(t) sum to a boundary value as before and only the 
weight corresponding to ki contributes in W{t). We will encounter such a situa- 
tion for driven lattice gases. 

3 Bulk driven lattice gases 

An important example of an SNS with spatial structure is a stochastic lattice 
gas where particles jump at random to neighboring lattice sites. We envisage 
two mechanisms for how the conserved particle density is driven away from equi- 
librium. The first one, studied in this section, is a global driving force which 
corresponds to driven diffusive systems (DDS). In the second one, considered in 
the next section, particles are injected/removed at the boundaries whereas the 
bulk is governed by reversible dynamics satisfying detailed balance with respect 
to an equilibrium stationary measure. In both cases the full dynamics does not 
satisfy detailed balance with respect to the stationary measure fj, s . 

We consider particles with exclusion hopping on a regular lattice 7h d . The 
configuration is specified by the occupation variables r\ x , x G Z d , taking the values 
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(empty) and 1 (occupied). Particles on different sites may interact by making 
the jump rate of a given particle depend on the configurations on nearby sites. As 
explained in Section 2, the fluctuation theorem holds in general. For notational 
simplicity we restrict ourselves however to the one-dimensional lattice, d — 1. 

We first study particles subject to a uniform drive. This will set up a steady 
state current in a finite system of particles jumping on a ring, x = 1,...,£, 
with x = I + 1 identified with x — 1. A particle configuration is denoted by 
7] = (rji, . . . , rji). We denote by c xx+ i(rj) the exchange rate for the pair of sites 
(x,x + 1), i.e., if T] x = 1 (0) and rj x+ i = (1), then c xx+ i{rf) is the rate at which 
the particle at x (x + 1) jumps to the empty site at x + 1 (x). 

In equilibrium the distribution of particles is determined by the Gibbs measure 
H eq ~ exp[— /3H] where (3 is the inverse temperature and Hirf) is the energy 
function. The rates c xx+ \ satisfy the condition of detailed balance, with respect 
to ji eq whenever 

e- w <Wi(l) = e-^ xx+ \ xx+1 ( V M ) . (3.1) 

Here r] xx+l denotes the configuration rj with the occupations at x, x + 1 inter- 
changed. A simple example of H(r]) is the energy of a system with nearest 
neighbor interactions 

i 

H(T]) = - J ^T] X T] X+1 , (3.2) 
x=l 

but such an explicit form is not used here. Rates satisfying (3.1) lead the system 
to equilibrium. In order to drive the lattice gas by a uniform external force field 
of strength F we add to H(rj) the linear term — ^2 X Fxrj x . While this term does 
not respect the boundary conditions, energy differences always do. Therefore the 
rates for a uniform drive are assumed to satisfy 

c F xx+ M = c F xx+ i(v xx+1 ) exp[-/?(tf (rT +1 ) - H( V ) - F( Vx - Vx+1 ))] . (3.3) 

We have called this condition local detailed balance [|TT|, [T3[], because locally a 
particle feels a linear potential. It is only through the periodic boundary condi- 
tion, that the steady state becomes a nonequilibrium one. For a closed system 
the particles would pile up at the right (for F > 0). 

To state the fluctuation theorem we merely have to apply the results of Section 
2. The logarithmic ratio in (2.5) is 

- (3(H( V XX+1 ) - H( V )) + PF( Vx - Vx+1 ) (3.4) 

with f] the configuration before and r] xx+l the configuration after the jump. The 
energy difference summed over the jumps is telescoping and yields a pure surface 
term as in (2.36). Therefore the relevant part of the action functional is given by 



10 



W{t)=(3F [ *Tj x>x+l (s)ds. (3.5) 

Jx,x+i( s ) is the actual particle current across the bond (x,x + 1), i.e. for a given 
history {r](s), < s < t}, J x , x+ i(s) is a sequence of (^-functions located at the 
times of jump through the bond (x, x + 1) with weight +1 (—1) if the jump is 
from x (x + 1) to x + 1 (x). The action functional is f3F times the the total 
current integrated over the time span t, which equals the signed number of all 
jumps up to time t. 

The entropy production can be read off from (2.30). If we set nirf) = 

g(r])Z~ 1 e~ l3H ^ and denote by {-) eq the average over Z~ l e~^ H ^\ then 

Riti) = ^(^(^[e-^-^)^^ 1 )-^^)] 



2 = 1 



[loge-^"+^(rT +1 ) - logg(r))]) eq . (3.6) 



We clearly have R(fJ>) > unless F = |28| . The identification of R(n) as entropy 
production is further supported by taking /i to be a state of local equilibrium. 
In the limit of slow density variations R coincides with the phenomenological 
entropy production based on the nonlinear diffusion equation as a macroscopic 
equation for the density p8 . 



We emphasize that the action functional is equal to the total current only if 
the condition (3.3) of local detailed balance holds. To illustrate this point we 
consider nearest and next nearest neighbor jumps. Local detailed balance for 
next nearest neighbor jumps means 

C F XX+ 2(V) = c F xx+ M x+2 ) eM~P(H(v xx+2 ) ~ H( V ) - 2F( Vx - Vx+2 ))} . (3.7) 

As before, W(t) constructed according to the rule (2.5) equals the total current. 
There are many ways to violate (3.7). Just as an example, if in (3.7) we replace IF 
by aF, then the action functional is W{t) = j3F J dt^2 x (J x x+ i(t) +aJ X;X+ 2(t)) 
which equals the particle current only if a = 2. Thus, while the fluctuation 
theorem holds in great generality, the specific physical interpretation is linked to 
local detailed balance. We could also in this case consider the generator as the 
sum of two generators, one for nearest and one for next nearest neighbor jumps, 
and have more general fluctuation theorems as in Section 2.3. The physical 
interpretation would however remain obscure. 



4 Boundary driven lattice gases 

We use the same set-up as in the previous section. Only now the box consisting of 
sites {1, 2, . . . , £} is not periodic, and there is no driving force. Instead particles 
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are injected and removed at the left (x = 1) and right (x = £) boundary. The 
generator, L, for the full dynamics is naturally decomposed as a sum of three 
pieces, 

L — Lb + Li+ L£. (4.1) 

Lb is the dynamics inside the box satisfying detailed balance with respect to e~^ H , 
cf. (3.1) with x = 1, ...,£ — 1. Lb conserves the number of particles. Li models 
the particle source/sink at the right and left boundaries and has the form 

LJ( V )=c l ( V )[f(v l )-f(v)}, i = l,l. (4.2) 

Cj (77) is the rate for the transition rji to 1 — rji while the remaining configuration 
is untouched. The configuration after this transition is denoted by if . As before 



we will require the rate Ci{rj) to satisfy local detailed balance in the form [O 



Cifa) = exv[-(3{H{rf) - H(r,)) + /3 Ti (l - 2 Vi )) , (4.3) 

where local refers now to the boundary points which are coupled to reservoirs 
with chemical potentials 71 and 7^. Thus if we consider the dynamics generated 
only by L& + Li, then the dynamics will satisfy detailed balance with respect 
to the stationary distribution Z~ l exp[— (3H(rj) + f3jiN(j])] (achieved in the long 
time limit) with N(t)) = Ylt=i Vx the number of particles. The same will be true 
for the full L in (4.1), if 71 = 7^. If on the other hand 71 7^ 7^, then the sources 
are unbalanced and there will be in the steady state a net flux of particles through 
the system going from the reservoir with high to the reservoir with low chemical 
potential, i.e. from left to right if 71 > j£. 

We now apply the results of Section 2 (with Aj = A). The logarithmic ratio 
(2.5) reads 

- (3(H(^ +1 ) - H( v )) -PY, W) - # fa)) + E W) - iV (^))' ( 4 - 4 ) 

i=l,£ i=l,£ 

where 77 is the configuration before and r] xx+1 , resp. 77*, the configuration after 
the jump. Again the energy differences sum to a surface term. Let Ji(t) be the 
boundary current, i = 1,£. It is a sequence of delta functions located at times 
when r\i jumps to 1 — 77^ with weight 1 (—1) for the transition from to 1 (from 
1 to 0). When the boundary current J\ is positive there is a net flux of particles 
from the ith reservoir into the system. Then 

W(t) = (3 [ ( Tl Ji(s) + i t Ji{s))ds . (4.5) 
J 

Because of particle conservation there are many equivalent expressions for 
W(t) in the sense that they differ only by surface terms and have therefore the 
same rate function. Let us denote by N[ XtV ](t), x < y, the number of particles in 
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the interval [x, y] at time t. By the conservation law for the number of particles 
we have t 

N lx , y] (t) - N [x , y] (0) = [ (J x - llX (s) - J y>y+1 (s))ds . (4.6) 

Jo 

Since N[ x>y ^(t) < \y — x\ + l, J* * J x -i,x( s )ds and J Q * J yjy+1 (s)ds must have the same 
large deviations. If we include in the mass balance also the boundary currents, 
then we can replace in (4.5) Ji(t) by any J x , x +i(t) an d Jg(t) by any — J x , x +i(t) 
without modifying the rate function e(w). While this gives many equivalent 
choices, the simplest one is perhaps 

W x (t) = /%i - 7,) / J x , x+1 (s)ds = /3( 7l - lt )W x {t) . (4.7) 
Jo 

W x (t) differs from W(t) only by a surface term. The action functional W x (t) is 
the time integrated particle current across the bond (x,x + 1). The generating 
function for W x (t) satisfies the fluctuation theorem in the form e(A) = e(/3(7i — 
ji) — A), with e independent of the choice of x. 

The conservation law (4.6) can be used also in the driven lattice gas of the 
previous section to produce equivalent action functionals. For example W(t) of 
(3.5) is equivalent to 

W(t)=£pF f J x , x+1 (s)ds. (4.8) 

However the symmetric form (3.5) seems to be more accessible to analytic com- 
putations, as can be seen from the examples described in the Appendix. 

5 Fluctuation theorem for diffusion processes 

We consider next a general diffusion process Xt € MP. It is specified by a drift 
vector c(x) (a vector field on R") and a positive definite diffusion matrix a(x) > 
(an n x n matrix valued function on M. n ). For later convenience we write the 
generator as 

L = iv-(aV) + c-V (5.1) 

with V representing differentiation with respect to x. x t is the solution of the 
stochastic differential equation 

dx t = (c + ^(V • a))(x t )dt + y/a(x t )db(t) (5.2) 

with db(t)/dt standard white noise. We use Ito's convention for stochastic differ- 
entials. 
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We claim that for diffusion processes the action functional W(t) is given by 
W[t)=2 [ (a^cixs)) ■ dx s + f aV{a- 1 c)(x s )ds . (5.3) 



We may think of W(t) as the limit of the symmetrized (Stratonovich) approxi- 
mation, 

t/e 

W{t) = lime^(a- 1 c(x (j+ i )e ) + o -1 c(a: ie )) • (x (j+1)e - x jt ) . (5.4) 
e ^ 3=0 

If a is diagonal and independent of x, then a#/2 corresponds physically to the 
temperature In this case W(t)/{3 equals the work done by the force c on the 
system during the time span t, in accordance with the physical meaning of the 
action functional for lattice gases. 

As before, we expect that the equilibrium forces do not show up in the large 
deviations of W(t). For diffusion processes detailed balance with respect to e~ u ^ x ' 
means that 

c=--aVU. (5.5) 



Inserting into (5.3), by Ito's lemma ||29|| , 

- / VU(x s ) -dx s -- [ aWU(x s )ds = U(x ) - U(x t ) , (5.6) 
Jo 2 Jo 

which is indeed a pure surface term. Only driving forces make a contribution to 
W(t) proportional to t. 

Since a > by assumption, we are in the uniformly elliptic case, where the 
transition probability, e Lt , has a density and (e Lt ) X)X > > for t > 0. We require 
that the drift c is sufficiently confining. Then xt has a unique stationary measure 
with density fi s (x) > 0, satisfying \i s e Li = fi s . 

Let us denote by (■) expectation with respect to the stationary process x t 
(starting in the invariant measure /i s ). Then 

lim --log(e- w «) = e(A) (5.7) 

t— >oo t 

and the fluctuation theorem 

e (A) = e(l - A) (5.8) 

should hold. 

To prove (5.8), as for jump processes, we first define the function g(x,t) = 
K x (exp[— \W (t)}) as the expectation with respect to the stationary diffusion pro- 
cess Xt conditioned to start at x. Then g(x,t) satisfies 

^g(x,t) = L x g(x,t) (5.9) 
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with 

L\ = -V ■ (a ■ V) + c • V — 2Ac ■ V — A(V ■ c) — 2A(1 - A)c • a _1 c , (5.10) 
where the last two terms act as multiplication operators. We conclude that 

(e" ww ) = J j v s (x)(e L >% x ,dxdx'. (5.11) 
On the other hand, from (5.10), 

L* x = Li- X , (5.12) 

which implies (5.8) by Perron- Frobenius. 

Eqs. (5.10) and (5.11) can also be viewed as an application of the Cameron- 
Martin-Girsanov formula |[29|| . The path measure of the process generated by L\ 
has a density relative to the path measure of the process generated by L. We 



write this density as exp[i?]. Then, according to [^9| , Section 6.4 



R = -2A J a~ l c{x t ) ■ dx t + 2A J {a~ l c) ■ (c + ^Va)(x t )dt 

-2A 2 y c ■ a' 1 c(x t )dt - 2X(1 - X) J c ■ a _1 c(x t )dt - A J V ■ c{x t )dt 
= — 2A y a _1 c(x 4 ) • <ixt — A J aS/{a~ l c){x t )dt 

= -XW(t), (5.13) 
where we used the identity 

a(a)— a _1 (a) = — ( — a(a) J a _1 (a) (5-14) 



da \da 

valid for a parameter-dependent matrix a(a). 

From the argument of Section 2, we conclude that W(t)/t has large deviations 
as stated in (2.11) with a rate function satisfying the symmetry (2.10). 

6 Particle under mechanical and thermal drive 

In the stochastic differential equations of Section 5 the diffusion acts everywhere 
on the configuration space M n . This is appropriate for strongly overdamped 
stochastic systems such as those considered in time-dependent Ginzburg-Landau 
theories. In many physical situations, however, stochasticity is assumed to act 
only in velocity space. Then the matrix a in (5.1) has zero eigenvalues and, 
the results of Section 5 are not directly applicable. To understand the required 
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modifications we consider, as an example, a mechanical particle of mass m with 
position x in the d- dimensional torus T d and velocity v in M. d subject to noise 
and friction. The particle has the mechanical energy H = \mv 2 + U(x), where 
U is some periodic potential. In equilibrium the stationary distribution is given 
by exp[—(3H(x,v)] with the temperature. We envision two mechanisms to 
drive the system out of equilibrium: (i) There is a mechanical driving force F(x). 
The standard example is a constant electric field E, i.e. F(x) = E. If F has 
a periodic potential part, it may be added to U but as we will see this has no 
effect on the large deviation of the action functional, (ii) There is a nonuniform 
temperature The case of a periodic variation in (5 (and constant friction) 

has been studied recently in the context of Brownian motors |3(| |3l| with the 
at first sight surprising result that the steady state maintains a nonzero current. 
The dynamics of the particle is governed by the stochastic differential equation 



d 2 

m dt^ Xt = ~^ u ( x t) + F ( x t) - mj(x t )v t + (2mj(x t )/P{x t ))^{t) . (6.1) 



Here 7(2:) > is the friction coefficient, (3{x) > the inverse temperature and 
£(t) standard white noise. The generator of the corresponding Fokker-Planck 
equation reads 

L = v ■ V x + \--VU + —F — 7(x)v J • V„ + (rf(x)/mP(x))Vl . (6.2) 



mm 

We claim that the action functional for the fluctuation theorem is given by 

W(t)= [ (3(x s )F(x s )-dx s + [ H(x s ,v s )df3(x s ) (6.3) 
JO Jo 

for x s ,v s a solution to (6.1). The first term is /3x(work done by F) and the 
second term represents f3 x (work due to thermal gradients) integrated along the 
trajectory x s , < s < t. If F — — VV, then up to a surface term 
Jq P(x s )F(x s ) ■ dx s = j^V{x s )df3{x s ) which, in (6.3), should be added to the 
energy H(x, v). 

Let (•) denote the stationary average for (6.1). Then 

lim --log(e- ww ) = e(A) (6.4) 

t— »oo t 

and — e(A) is the maximal eigenvalue of 

L x = L - \(p{x)F{x) ■ v + H(x, v)v ■ V x (3(x)) , (6.5) 

the last two terms being considered as multiplication operators. The correspond- 
ing maximal eigenvector is denoted by f\(x, v) > 0. Let R be the velocity reversal 
operator, Rf(x,v) = f(x, —v). Then 

Re~^ x)H L x e^ x)H R^Re-^ 11 f x (x, v) = -e(X)Re~^ x)H f x (x, v) . (6.6) 
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By a straightforward computation 



Re-^L.e^R-^Ll^. (6.7) 

Therefore 

L\_ x Re~ pH h = -e{\)Re- pH h = -e(l - X)Re^ H f x (6.8) 

by Perron- Frobenius, since Re~^ H f\(x, v) > 0. We conclude that the fluctua- 
tion theorem holds and e(A) = e(l — A). 

As for lattice gases, we expect that the action functional W(t) is linked to 
the entropy balance. To spell out the details we consider the time-dependent 
probability density p t (x,v) of (x t ,v t ). It satisfies J^p t = L*p t . As usual the 
Gibbs entropy is given by Sc(pt) = — f Pt log ptdxdv and changes in time as 



j t S G ( P t) = °(pt) + Upt) ■ (6.9) 



The first term 



^—[>/mf3vpt + —=V v p t ] 2 dxdv (6.10) 
Pt Vm(3 

is positive definite. We identify it as the entropy production in the system. The 
remainder reads 

js(pt) = J PtPHdxdv — J ptf3F ■ vdxdv — J ptHv ■ V x f3dxdv . (6.11) 

We regard it as the entropy flow from the system to the mechanical and heat 
reservoirs. 

According to the rules used before, we construct W(t) by integrating the 
linear functional defining the entropy flow along a stochastic trajectory. Using 
(6.11) this yields 

~PH{x„ v s ) + pF(x s ) ■ v s + H{x s , v s )v s ■ VP{x s ) J ds . (6.12) 
Since the first summand is a surface term, we recover (6.3). 



7 Stochastic and thermostatting heat reservoirs 

Physical systems with stationary flows of heat, momentum, etc., are usually 
modeled by coupling the system to thermal reservoirs represented by stochastic 
forces acting near the boundaries [0-0. Away from the boundary the time 
evolution obeys the same laws as an isolated system. To be specific, let us 
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consider one such model of heat flow through a classical system of iV particles 
of mass m with position qj, momentum pj,j = 1, . . ., N, inside the slab A = 
[-£-a, £+a] x [0, if. Let the boxes A_ = [-£-a, -i] x [0, if, A+ = [£, i+a] x [0, if 
be the left and right boundary zones. The indicator functions of these sets are 
denoted by \±. The particles interact through a short range pair potential V and 
are confined to the slab by the wall potential V w . The equations of motion are 



d 1 

if = m Pi > 

" Pj = F j + F w (q j ) + J2x S (q j )(-lP j + (2m 1 /(3s) 1 %(t)), (7.1) 



dt 



8=± 



where qj G A,pj G M. 3 , Fj is the force acting on the j-th particle, Fj = 
— / , =1 W(qj — qi),F w = —W w is the force from the wall. are a 

collection of independent white noises. The friction and the stochastic forces 
operate only when the particle is in the boundary layers acting there like thermal 
reservoirs at inverse temperature (3$. If P+ — P — then the steady state is 
Z^expl-PH), H = + VM) + |E£=i - Qi)- If 0+ * 

there is, in the steady state, a constant heat flux from the higher to the lower 
temperature reservoir. (This is a variation of the model used to describe heat 
flow in a crystal where the Langevin forces act on the particles in the end layers 

0- 0)- 

We now follow the method explained in Section 6. For the similarity trans- 
formation we use 



exp 



N f 1 1 N 

- E ft*) \ k$ + v «(<to + 2 E v (* " <&) 

3=1 I i=l 



(7.2) 



The local inverse temperature (3{q) = (3s for q G As and /3(q) interpolates 
smoothly otherwise. For each choice of (3(q) we obtain a distinct action functional 
W(t). By local conservation of energy, these differ only through boundary terms. 
For later purposes it is convenient to choose a step interpolation as (3{q) = (3- 
for qi < and (3(q) = (3 + for q\ > 0. Then the action functional is 

W(t) = {(3 + - /?_) / ds J (s) (7.3) 

with J the energy current through the plane {gi = 0}, 



N 1 f 1 1 * 1 

J o = E -a- 2^ + 2 E ^ - «o 

i^i r 1 

i 1.7 — 1 



(7.4) 
}1 + (1 - A)%i) , 
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where q^ is the x-component of the jth particle position vector. 

As an advantageous numerical alternative deterministic thermostatting reser- 
voirs have been developed [[RJ. We mention here the proposal of Gallavot-ti and 

The equations of motions are Newtonian in 
and A_ additional friction terms are added 



Cohen [22] for thermal reservoirs, 
the bulk. At the boundary layers A_| 
which are constructed in such a way that the total kinetic energy of the particles 
in A + , resp. A_, is kept constant. The particular boundary temperature is then 
fixed through the initial conditions. The equations of motion read 



d_ 

It 



1 



m 



Pj 



dt 



Pj 



Fj + F w (q j ) + Y^ Xs(qj)asPj 



(7.5) 



s=± 



The 'friction' coefficients are given by 



j=i 

+- X s(q j )(F j + F w ) ■ Pj \ , 6 = ±1. 



(7.6) 



For thermostatted systems the action functional is the phase space contraction 
integrated along a trajectory of the dynamical system, which in our case equals 



3Na = 3N + a + + 3ALa_ + 0(1) 



(7.7) 



with N± the number of particles in A±. The fluctuation theorem holds provided 
the dynamics is sufficiently hyperbolic for the SNS to be described by an SRB 
measure. This is not expected to be true in general. It is then an assumption, 
embodied in the chaotic hypothesis [EM], that the fluctuation theorem remain 



valid for realistic physical systems. The numerical simulation in j32] gives strong 
support for the validity of the fluctuation theorem in one such model. There the 
two heat reservoirs are linked through an anharmonic chain, rather than a fluid. 
Further numerical support comes from a study of shear flow for hard disks in a 
box [Q where momentum but no energy is transferred to the system at opposite 
boundaries. The fluctuation theorem is well verified 0], although it seems to be 
difficult to get beyond the quadratic approximation. 

The action functional defined with the phase space contraction (7.7) appears 
to be rather different than the one of (7.3). However the conservation law can be 
used again to transform a boundary flux to an interior flux. Let Hi be the total 
energy of particles in the left half box [—a — £, 0] x [0, £} 2 . We have 



d_ 

7ft 



N 



Hi{t) 



-Mt) - 



(Qj)p 2 j 



a. 
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since Ysj=iX-{Qj)Pj * s a constant of motion and initially fixed to 3iV_//3_. We 
conclude that 

-t rt 
ds3Na(s) = -H(t) + H(0) + (/?+ - / dsJ (s) . (7.9) 

Jo 

Up to a boundary term the phase space contraction is just the energy flux across 
the plane {q± = 0}, the same quantity which has been obtained for stochastic 
reservoirs, compare with (4.7). 

It is tempting to assume that for stochastic and thermostatted boundaries, 
the rate functions of the fluctuation theorem are identical. This requires however 
a very strong form of equivalence of ensembles, in the sense that in the interior 
of the system even the large deviations of the current are in the steady state 
independent of the mechanism by which the boundaries are cooled/heated. This 
appears to be the case for the model studied in p2| . On the other hand recent 



simulations |36) suggest that this equivalence fails for the model described in 



8 Green-Kubo formula, Onsager reciprocity 

As observed in [24| the fluctuation theorem yields the Green-Kubo formula at 



equilibrium. We give an alternate derivation in the spirit of po| . For the sake 
of concreteness we consider a fc-species lattice gas on a <i-dimensional lattice. 
The lattice gas is driven by an external field, cf. Section 3, and satisfies local 
detailed balance. For each species the driving field is a <i-component vector. 
It is convenient to regard the driving field as the m-component vector F = 
(Fi, . . . , F m ), m = d x k. There are m different types of jump, labeled by the 
species and lattice directions. Therefore the generator of the lattice gas is the 
sum 



5>w 



3=1 

where each L^> satisfies the conditions of Section 2.3, compare with (2.23). To 
each there is associated the current Jj(t). 
We define 

m rt 

e(F] Ai, . . . , A m ) 



lim--log(exp[-/3V / A,- J^s) ds])p , (8.2) 



where we have, in the action, absorbed Fj into the X/s and indicated explicitly 
the dependence of the stationary state and the dynamics on F. The fluctuation 
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theorem now takes the form 



e(F*; Ai, . . . , A m ) = e(F*; F t - A 1; . . . , F m - A m ) . (8.3) 

We have (Jj(t))o = 0, since F = corresponds to equilibrium. Differentiating 
(8.2) at A = 0, we obtain 

() e{F*;0)=P(J j {t)) p , (8.4) 



the average j-th current. 

Now the linear response in the average current to a small driving field is given 

by 

w^ eM = aii - (8 ' 5) 

Differentiating (8.3) with respect to F and A at F = A = we conclude that 

/3aij = -f3<Jij + f3 2 &ij (8.6) 

and therefore 

°ij = -Ptij , (8.7) 

where 

This is the standard Einstein-Green-Kubo relation between linear response in the 
current and the time-integrated current-current correlation, which satisfies the 
Onsager relation = er^. 

We remark that the derivation of (8.7) differs considerably from the standard 
computation, we refer to [28] for an exposition of the latter. There one differ- 
entiates, at finite volume, the steady state current (jp) p at F = 0, which by 
definition equals er^-. This gives two terms: one from differentiating the current 
function jg and one from differentiating the steady state (•) g. Their sum is then 
(8.8), because ( Ji(t) Jj(0))q has a 5-peak at t — 0, whose weight is the first term, 
and a smooth piece, whose time integral equals the second term. 



9 Conclusions 

Within the framework of stochastic dynamics the fluctuation theorem holds in 
great generality. For finite state spaces, like stochastic lattice gases, we have 
given a proof. For the nondegenerate diffusion processes of Section 5 we could 
provide sufficient conditions on the driving forces and the diffusion matrix which 
ensure the validity of (5.8). The situation is more delicate for degenerate diffu- 
sion processes, in particular for deterministic bulk systems driven by stochastic 
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boundaries. Since the noise vanishes on large parts of the phase space, even the 
convergence to the steady state is not obvious. To our knowledge there are only 
a few rigorous investigations of SNS of this type || f|, |, [7], [H| . In 0, classical 
point particles are considered. At a collision with the wall they are reflected 
with a Maxwellian distribution at the local temperature of the wall. Existence 
and uniqueness of the invariant measure is established provided the forces are 
repulsive and have a range of the size of the box. This latter property prevents 
pockets in phase space which never see the wall. In || explicit SNS measures 
are obtained for a harmonic chain coupled to stochastic reservoirs specified by 
Langevin forces at the ends of the chain. While in |10[ an anharmonic chain 



is coupled to a free field leading to similar boundary conditions. For this sys- 
tem existence and uniqueness of the invariant measure is established in under 
a Hormander condition. We expect that for these systems the GC fluctuation 
theorem can be proved. Thermostatted boundaries are beyond our mathematical 
abilities at the present time. 

Viewed from the physics point of view the fluctuation theorem is a conse- 
quence of time-reversal. For a given stationary stochastic process one considers 
the density of the path measure of the time-reversed process relative to the origi- 
nal one. The logarithm of this density is the action functional which satisfies the 
fluctuation theorem. At this level of generality it is somewhat unexpected that in 
concrete models, which satisfy the condition of local detailed balance, the action 
functional has a direct physical interpretation. By separating the rate of change of 
Gibbs entropy into a production and flow term we identify the observable whose 
average defines the entropy flow. The action functional is then this observable 
integrated along a stochastic trajectory. Generically it is the current, multiplied 
by the driving force, corresponding to the conserved field. The GC fluctuation 
theorem is then a symmetry property for the large deviations in the current. Such 
large deviations are not readily observed in physical systems, because they refer 
to exponentially small probabilities. Also in numerical experiments special efforts 
are needed [E3J, pi| , P5fl . Even worse, the action functional is extensive and the 
probabilities are exponentially small also in N. *Thus it would be desirable to 
have a local quantity satisfying (at least approximately) the fluctuation theorem 
but this does not appear to be the case for the systems investigated numerically 
thus far ||, @. 



There is one prediction of the fluctuation theorem which is worthwhile to 
emphasize. We use the convention of the previous section where the symmetry 
reads e(A) = e(F — A), i.e. e is even relative to A = F/2, and we may expand as 
e(A) = eo + |e2A(F — A) + . . .. Assuming that in the interval — 5 < A < F + 5 with 
some small 5 > the quadratic approximation is valid, then by differentiating at 
A = we obtain (j) F = e 2 F/2, where e 2 = JZo dt((J(t)J(0)) F - (J(t)) F (J(0)) F ) 
will generally depend on F. The nonlinear response in the current is then linked 
to the time-integral over the current-current correlation. The knowledge of e(A) 
thus provides information about the range of validity of the Einstein relation and 
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linear response theory. 



A Appendix: Large deviation function for the 
asymmetric exclusion process 

To determine the large deviation function in any of the models we have discussed 
is like computing a free energy in equilibrium statistical mechanics. Therefore 
it is not surprising that even in the stochastic framework explicit results are 
available only for a few models. We discuss here the large deviation function for 
the asymmetric simple exclusion process in one dimension where particles random 
walk with an average drift and interact through the constraint that there can be at 
most one particle per lattice site. As a warm up we consider this model without 
constraint, i.e. independent biased random walks on a lattice. It would be of 
interest to have also an example for the boundary driven lattice gases of Section 
4. Besides independent particles, the simplest model is symmetric exclusion in 
the bulk and boundary reservoirs as discussed in Section 4. Its large deviaiton 
function, for W x (t), satisfies then e(A) = e(/3(7i — 7^) — A), cf. (4.7). Thus in 
the limit 71 — > 00, i.e. at the left boundary particles are only injected, the GC 
fluctuation theorem predicts lim t ^ 00 (— 1/t) logProbdW 7 ^ < at}) = for a < 0. 

We first consider N independent particles, i.e. no exclusion, on a ring of size 
I driven by a uniform force of strength F . Since particles are independent, it 
suffices to study the one-particle problem. For N particles we merely have to 
multiply by N at the end. It is convenient to set (3F/2 = E. The generator for 
the dynamics reads 

L E f(x) = (e E /2 coshE)f(x + 1) + (e~ E /2 coshE)f(x - 1) - f(x) . (A.l) 

The elementary jump time is normalized such that for E — > 00 the jump rate 
to the right is one. W(t) is now the total number of signed jumps up to time t. 
To compute the generating function (e~ xw ^) for large t we need the maximal 
eigenvalue of 

L E ,xf(x) = (e E - x /2 cosh E)f(x + 1) + ( e - E+x /2 cosh E)f(x - 1) - f(x) . (A.2) 

Clearly / = 1 solves L F \f = —e(X)f with 

e(A) = 1 - (cosh(£ - A)/ cosh E) . (A.3) 

Its Legendre transform is given by 

e(w) = —Ew + 1 + wargsinh(w cosh E) — a/w 2 + {coshE)~ 2 . (A. 4) 

e vanishes at the average current j = tanh E. On general grounds we know that 
e(A) is convex down, e(0) = 0, e'(0) = — j, the average current, e(A) = e(2E — A) 
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and that e(w) is convex up, e(w) > , e(wo) = for wo = j = j(W(t)), 
e(w) — e(—w) = —2Ew. 

Note that in the limit E — > oo, e degenerates to e(w) = 1 — w + wlogw 
for w > and e(u>) = oo for w < 0, which merely reflects that left jumps have 
probability zero for E = oo. 

We expect more structure for particles interacting through hard core exclu- 
sion. This corresponds to H{rj) = in the notation of Section 3. The jump rates 
with a uniform driving field are then 

= O^/ 2 cosh E )Vx(l -Vx+i) + (e- E /2 C oshE)(l-r] x )r] x+1 . (A.5) 

We normalized the rates such that for E — > oo a particle jumps with unit rate to 
the right. (A.5) are the jump rates of the asymmetric simple exclusion process 
(ASEP). 

The action functional W(t) is the total current integrated over the time span 
t. Expectations, (•), are taken in the stationary process at a given number of 
particles N, 1 < N < £ — 1. For the ASEP this means that at time t = all 
allowed configurations have equal weight. Then 

lim -- log(e- w W) = e E (X, £, N) (A.6) 

t^oo t 

and — e E is the maximal eigenvalue of L\ with 

e 

L x f( V ) = (2cosh E)- 1 {(e E -\(l - Vx+i) + e^ E - x \l - Vx ) Vx+1 ) 

x=l 

f(v xx+1 ) - (e E Vx(l ~ ife+i) + e~ E (l - Vx)Vx + i)m} ■ (A.7) 

The fluctuation theorem reads 

e E (X) = e E (2E - A) , (A.8) 

which can be seen also directly from (A.7). 

L\ is a linear operator with a structure which is similar to that of the hamil- 
tonian for a quantum spin chain. To make contact with that model we rewrite L\ 
in terms of the Pauli spin \ matrices a x = (o - ^, a x1 a x ) associated to every lattice 
site x. We identify i] x = 1 with spin up, a\ — 1, and i] x = with spin down, 
a\ = —1. Then 

i 

-L X = H X = J2 {(V4)(l - °y x+1 ) - (e E - X /2 cosh Ejfc 

x=l 

-(e- E+x /2 cosh E)aZat +1 ) , (A.9) 

where = |(<7^±icr^) are the spin raising and lowering operators. We note that 
HI 7^ H\, unless E = A, and H\ does not have directly a quantum mechanical 
interpretation, except through analytic continuation in A. 
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To determine the ground state energy of if a, the only technique available is the 
Bethe ansatz. Fortunately, H\ is covered by the famous paper of Sutherland and 
Yang [57]. However, the analysis of the nonlinear Bethe equations still require a 



considerable effort, which in our section of parameter space has been carried out 
only recently. We use the results by D.Kim [Q. He is in fact interested in the 
energy gap, but also gives e E (X, £, N) to leading order in N in the limit I — ► oo, 
N — > oo, p = N/£ fixed. We expect that also order 1 terms could be obtained, 
but such an analysis is certainly beyond the scope of our paper. In the limiting 
case E — > oo, the Bethe equations simplify. This has been exploited recently by 
Derrida et al. [^, |40| . They determine eoo(A,£, N) fairly explicitly and discuss 



the order 1 corrections to the infinite volume limit. 
To leading order in N, N/ i = p, we have 

e E (XJ,N)=£e E (X,p), (A.10) 

i.e. e E is extensive. e E has a flat piece, e E (X, p) = for < A < 2E. At A = 0, e E 



has the slope p(l — p) tauhE and the asymptotics |38[ 

e E {X,p) = p(l -p)(tanh£)A (A.ll) 

(4p(l - p)) 4/3 (tanh£)|A| 5/3 (l + C(|A| 2/3 )) 



1 /Q \ 2 / 3 
1 / 07T 
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for A < 0. By symmetry e E has the same asymptotics at A = IE. For E — > oo, 
the right zero e^(A, £, AT) moves to oo and the flat piece of e E (X, £, N) never levels 



off, in agreement with |39[ 40 



e'(0) is just the average current which is p(l — p)ta.nhE. e"(0) is formally 
the integral over the total current- current correlation and therefore the mobility 
of the lattice gas [p8 |. In our case e"(0) = oo, which reflects the fact that in the 



ASEP density fluctuations propagate superdiffusively [fil 

At E = (A. 6) gives the large deviations in the total current at equilibrium, 
i.e. for the symmetric simple exclusion process. The fluctuation theorem reduces 
then to eo/,Ar(A) = eo,£,jv(— A), which also follows from time reversibility. eo(A,p) 
has no longer a flat piece. For small A, eo(A,p) = — |A 2 p(l — p). Thus the 
mobility a(p) = p(l — p) and the diffusion constant, which is o divided by the 
compressibility, D(p) = 1, independent of the density in agreement with the 
known bulk diffusion of the symmetric simple exclusion process. The more precise 



asymptotics |38|] is given by 



eo(A,p) = -ip(l-p)A 2 -l2^(2^/3( p (i_ p ))V3| A |8/3 (1 + 0(|A| 4/3 ) , (A12) 

On a formal level the "Burnett coefficients" correspond to higher derivatives of 
e E (X) at A = ||42|| . (A. 12) shows that for the symmetric simple exclusion process 
the Burnett coefficients are infinite. 
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Just as free energies, the large deviation function may become singular in the 
infinite volume limit and this happens already for the first nontrivial example. 
To understand the origin of this behavior, it is necessary to analyze those con- 
figurations which give the main contribution to (e~ xw ^). If A — ► oo only left 
jumps are permitted and similarly for A — > — oo only right jumps are permitted, 
essentially independent of the driving field E, E > for the sake of discussion. 
However, in the intermediate regime < A < 2E, the particles have two op- 
posing instructions. Should they follow E or A? In fact, neither. As can be 
inferrred from the maximal eigenvector, they just stick and form one big cluster 
43|| . The clustering is seen most easily at A = E, where H\ is the ferromagnetic, 



anisotropic Heisenberg model and the maximal eigenvector its ground state. At 
our parameters, H\ is in the ferromagnetic phase. Since the magnetization is 
fixed, typical configurations condense into one large cluster. In our units the 
corresponding ground state energy is of order 1 independent of £. 

It is instructive to reconsider our result from the point of view of the prob- 
ability distribution of W(t)/t, which is obtained from the Legendre transform 
of eg(A). To make it well defined, we note that, for large N, e^(A) smoothly 
interpolates between e#(0) = = e^(2_E), such that e' E (0) = Nj(p),e' E (2E) = 
—Nj(p),j(p) = p(l — p)tanhE, and ce{E) = 0(1). Thereby we obtain, for 
E > 0, 

eE ^ p) = \-2Ew for -j(p) <w<0. (A ' 13) 

For w slightly larger than j(p) we have e^(w,p) ~ (w — j(p)) 5 ^ 2 and eE(w,p) 
~ w log w in the limit of large w. For w < —j{p) the values of ce are fixed by 
the fluctuation theorem. 

To understand the iV-dependence of 6e(w, £, N)/£, we note that for W(t)/ 
Nt to have a value larger than j(p) we have to speed up all iV particles. For it 
to have a value smaller than j(p) it suffices to slow down a single particle, since 
the other particles pile up behind. Typical configurations consist of essentially a 
single cluster. The large deviation rate is 0(1/N) which on our scale corresponds 
to e-E = 0. Such a mechanism works only for (W(t)/Ni) > 0. To have a negative 
total current, again all iV particles are forced to move, now opposite to E. The 
linear decrease of 6e{w < N,£)/£ with corrections of 0(1/N) for —j(p) < w < 
follows from the fluctuation theorem. In the limit E — + oo, ce{w,£,N) = oo for 
w < 0. This just reflects the fact that the underlying process does not allow for 
histories with negative currents. 

From our discussion of e we conclude that the non-smooth behavior of e(A) for 
large N comes from blocking through slow particles, which suggests that, if we 
go to higher spatial dimension or soften the hard core, a smooth large deviation 
function may be recovered. But this remains to be seen. The symmetric case, 
E — 0, warns against too fast a conjecture. 
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